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Ab-initio study of BaTiC^ surfaces 

J. Padilla and David Vandcrbilt 
Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08855-0849 

(February 18, 1997) 

We have carried out first-principles total-energy calculations of (001) surfaces of the tetragonal and 
cubic phases of BaTiOs. Both BaO-terminated (type I) and Ti02-terminated (type II) surfaces are 
considered, and the atomic configurations have been fully relaxed. We found no deep-gap surface 
states for any of the surfaces, in agreement with previous theoretical studies. However, the gap is 
reduced for the type-II surface, especially in the cubic phase. The surface relaxation energies are 
found to be substantial, i.e., many times larger than the bulk ferroelectric well depth. Nevertheless, 
the influence of the surface upon the ferroelectric order parameter is modest; we find only a small 
enhancement of the ferroelectricity near the surface. 
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I. INTRODUCTION 

Recently there has been a surge of interest in the 
application of first-principles density-functional calcula- 
tions to the study of the rich phenomenology of the per- 
ovskite oxides, with special attention to ferroelectric (FE) 
properties!!! From these investigations, it has been found 
that the FE instability in these materials occurs as a re- 
sult of a delicate balance between long-range Coulomb in- 
teractions that favor the FE state, and sJaprt-range forces 
that favor the cubic perovskite phase.BO Moreover, the 
ferroelectric properties are well know to degrade in thin- 
filma and particulate^ geometries, suggesting that the FE 
state could be very sensitive to surface effects. 

The cubic perovskites have the chemical formula 
AB0 3 . For II-IV perovskites (e.g., BaTi0 3 ) A is a di- 
valent cation and B is a tetravalent transition metal, 
while for I-V perovskites (e.g., KNbOs) they are mono- 
and pentavalent respectively. The (001) and (111) sur- 
faces of these -materials have been most investigated 
experimentally^ There are two possible terminations of 
the (001) surface: the AO-terminated surface (type-I) 
and the B02-terminated surface (type-II). In II-IV per- 
ovskites, the AO and BO2 layers are charge-neutral, so 
that both type-I and type-II surfaces are non-polar. For 
I-V perovskites, the corresponding surfaces are instead 
polar. As for (111) surfaces, the atomic planes in this 
direction are of the form AO3 and B, and are charged in 
either case, so that the (111) surfaces are polar. Since 
polar surfaces are expected to be relatively unstable, we 
have chosen to focus here on the (001) surfaces of a II-IV 
perovskite, BaTiOs. 

Due in. part to the catalytic properties of SrTi03 and 
BaTi03,LJ there has been a continuous interest in the sur- 
face properties of these materials. There have been previ- 
ous theoretical studies especially for the case of the para- 
electric SrTiCw surface, but also for BaTiOs. Wolfram 
and coworkersQ, using a linear combination of atomic or- 
bitals (LCAO) cluster method, predicted mid-gap sur- 
face states for SrTi03, in disagreement with experimental 



investigations .BE] Only after ad hoc modifications to this, 
model could the. experimental results be accounted for .Ell 
Tsukada et alHl employed the DV Xa cluster method to 
study SrTiOs surfaces, finding no mid-gap surface states. 
However, cluster methods are not very suitable for high- 
accuracy calculations of relaxations and electronic states 
on infinite surfaces, underlining the need for the applica- 
tion of more accurate, self-consistent techniques. While 
such techniques have recently, yielded a great deal of in- 
sight into bulk perovskites ,crE3il3 their application to the 
study of perovskite surfaces has not been very extensive. 
In fact, we only know of two such studies. CoherEMZl 
presented linearized augmented plane wave (LAPW) cal- 
culations performed for slabs of tetragonal BaTiOs with 
(001) and (111) surfaces, using both symmetrical and 
asymmetrical terminations. Although some relaxations 
were allowed, the atomic positions were not fully re- 
laxed. Kimura— ft alxB used a plane-wave ultrasoft- 
pseudopotentiaO approach (as in the present work) to 
study the Ti02-terminated (001) surface of SrTi03, with 
and without oxygen vacancies at the surface. Again, the 
slabs were not fully relaxed. 

In contrast, we study here symmetrically-terminated 
type-I and type-II surfaces of tetragonal and cubic 
BaTiOs (001) for which the coordinates have been fully 
relaxed by minimizing the total energy. This allows us 
to study the influence of surface relaxation effects upon 
the FE distortion. For the tetragonal phase, we con- 
sider only the case of the tetragonal c axis (i.e., polar- 
ization) parallel to the surface. (Polarization normaJ to 
the surface is suppressed by depolarization fields J3 at 
least for clean surfaces.,)-. We employed the ultrasoft- 
pseudopotential methodta within the local-density ap- 
proximation (LDA). This technique permits us to calcu- 
late the Hellmann-Feynman forces on each atom, mak- 
ing it possible to find the relaxed structure much more 
efficiently than for methods that compute only total en- 
ergies. 

Experimental studies of perovskites surfaces are com- 
plicated by the presence of surface defects^ making it 



difficult to verify the surface stoichiometry. Therefore, 
most experimental investigations have not been very con- 
clusive. On SrTiC>3 surfaces, the situation is better: the 
surface relaxations have been studied£il and (as men- 
tioned above) the absence of mid-gap surface state has 
been demonstrated.0 For BaTi03 surfaces, the experi- 
mental reports_seem less conclusive. For example, evi- 
dence both fortJ and againsto surface gap states in this 
material have been published. 

Regarding the degradation of FE propertieSpfjor thin 
films and small particles as mentioned abovepB there 
does not seem to be any consensus about the origin of 
these effects. One possibility is that it is completely in- 
trinsic, i.e., that the very presence of the surface sup- 
presses the FE order in the vicinity of the surface. How- 
ever, there are many other possible causes. These include 
the effects of surface-induced strain; perturbations of the 
chemical composition near the surface related to the pres- 
ence of impurities, oxygen vacancies, or other defects; 
and the depolarization fields for the case of particles (or 
for films with polarization perpendicular to the surface) . 
Here, we take a modest step in the direction of sorting 
out these effects by characterizing the purely intrinsic 
coupling between the presence of a surface and the FE 
order, for the case of a free (vacuum-terminated) surface. 
As we shall see, we find the surface relaxation energies 
are large compared to FE distortion energies. However, 
we find very little effect for type-I surfaces, and only a 
modest enhancement of the FE order at type-II surfaces, 
with indications that it will be mainly confined to just the 
first few atomic layers near the surface. Thus, it appears 
unlikely that intrinsic surface effects are responsible for 
the degradation of FE order in thin films and particles. 

The remainder of the paper is divided as follows. In 
Sec. II, we describe the technical aspects of our first- 
principles calculations. In Sec. Ill, we report the results 
of our simulations. Finally, the main conclusions of the 
paper are reviewed in Sec. IV. 



II. THEORETICAL DETAILS 

We carried out self-consistent total-energy pseudopo- 
tential calculations in which the electronic wave func- 
tions were expanded in a plane-wave basis. The core 
electrons were frozen, and for a given geometry of the 
ions, the valence electron wavefunctions were obtained 
by minimizing the Kohn-Sham total-ienergy functional 
using a conjugate-gradient techniquerJ The exchange- 
correlation potential was treated with the LDA approxi- 
mation in the Ceperley- Alder form£3 The forces on each 
ion were relaxed to less than 0.02 eV/A using a modified 
Broyden scheme r3 ._. 

The Vanderbilt ultrasoft pseudopotential schemeO was 
employed. In this approach, the the norm-conservation 
constraint is relaxed, allowing one to treat rather local- 
ized orbitals with a modest plane-wave cutoff. The pseu- 



TABLE I. Computed and experimental values of structural 
parameters for BaTiOs in bulk cubic and tetragonal phases. 
a and c are lattice constants; <5 X are displacements associated 
with the FE instability as a fraction of c. Oi is the oxygen 
lying along x from the Ti atom, and <5 x (On) = <5 x (Oin) by 
symmetry) . 



Phase 


Parameter 


Exper. a 


Theory 10 


cubic 


a (A) 


3.996 


3.948 


tetrag. 


a (A) 


3.992 


3.938 




c(A) 


4.036 


3.993 




5 x (Ti) 


0.0135 


0.0128 




<5x(Oi) 


-0.0150 


-0.0150 




5x(Om) 


-0.0240 


-0.0232 




dopotentials for Ti, Ba and O are identical to those used 
previouslycJ in a study of bulk perovskites. The semi- 
core Ti 3s and 3p states and Ba 5s and 5p states are 
included as valence levels. A plane-wave cutoff of 25 Ry 
has been used throughout; previous work has shown that 
the results are well- converged at this cutoff. 

BaTiOs undergoes a series of phase transitions as 
the temperature is reduced, from the high-symmetry 
paraelectric cubic phase to FE phases with tetrago- 
nal, orthorhombic, and rhombohedral unit cells. The 
tetragonal structure is of special interest, since it is 
the room-temperature structure. In this paper, we are 
thus primarily interested in the surfaces of the room- 
temperature tetragonal phase, although for compari- 
son we also present results for surface of the elevated- 
temperature cubic phase. Ideally one would like to do 
this by carrying out ab-initio molecular-dynamics sim- 
ulations at the temperatures at which these phases are 
stable, but unfortunately this is not practical. Instead, 
we carry out ground-state (T = 0) calculations, but sub- 
ject to the imposition of the appropriate (tetragonal or 
cubic) symmetry in order to prevent the system from 
adopting the true rhombohedral T = structure. This is 
clearly an approximation, but we think it is a reasonable 
one. The computed ground-state structural parameters 
for the cubic and tetragonal bulk phases are given in Ta- 
ble Ml together with experimental values for comparison. 
(Again, the theoretical values are for T = structures 
with the appropriate symmetry imposed.) The computed 
lattice constants a and c for the cubic and tetragonal 
phases are 1 — 2% smaller than the experimental ones; 
this underestimation is typical of LDA calculations. We 
use the theoretical unit cell parameters in all calculations 
presented here. 

As shown schematically in Fig. 1, the periodic slab 
corresponding to the type-I (BaO terminated) surface 
contains 17 atoms (4 BaO layers and 3 Ti02 layers). 
Similarly, the type-II (Ti02 terminated) slab contains 18 
atoms (4 Ti02 layers and 3 BaO layers). For both cases, 
the slabs were thus three lattice constants thick; the vac- 
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FIG. 1. Schematic arrangement of layers in the 

BaO-terminated (slab I) and Ti02-terminated (slab II) su- 
percell geometries. Layers 1 and 7 are surface layers. 



uum region was two lattice constants thick. The z-axis is 
taken as normal to the surface, and the M z mirror sym- 
metry with respect to the central layer was imposed in 
all cases. For surfaces of the cubic phase, mirror symme- 
tries M x and M y were also preserved. For the tetragonal 
case, the polarization vector and thus the tetragonal c- 
axis were chosen to lie along x (parallel to the surface); 
here, M y symmetry was respected but m x symmetry was 
allowed to be broken. As mentioned earlier, the choice 
of polarization parallel to the surface is motivated by the 
fact that no charge accumulation results at the surface 
for this case, and thus no depolarization fields appear.tl 
The case of the tetragonal c-axis lying-perpendicular to 
the surface was considered by Cohen.E3iiIl 

The calculations were done using a (4,4,2) Monkhorst- 
Pack meshEO This corresponds to three and four k-points 
in the irreducible Brillouin zone (BZ) for the cubic and 
tetragonal supercells, respectively. To test the conver- 
gence with respect to k-point sampling, we repeated the 
calculation for a (6,6,2) mesh, finding that the surface 
energy differed by less than 3%. When the vacuum re- 
gion was enlarged to three layers in thickness, the surface 
energy changed by less than 4%. 

In order to study the relative stability of the two kinds 
of surface terminations, it is necessary to introduced ap- 
propriate chemical potentialsHj To simplify our analysis, 
we think of the independent constituents of the slab as 
being BaO and TiC>2 units. We define Ef to be the for- 
mation energy needed to make bulk BaTiC>3 from BaO 
and Ti02, per formula unit. Thus we have 



Ee — E] 



BaTiOa 



E T 



aO 



E 



TiQ 2 



(1) 



(by convention, E f > 0). E Bs ,Tio 3 , ^BaO, and E T io 2 are 
the energies of the bulk crystals, per formula unit, mea- 
sured relative to isolated ion cores and electrons. BaTi03 
was calculated in the relaxed tetragonal structure, and 



TABLE II. Atomic relaxations of the Ba-terminated sur- 
face (slab I) in the cubic (C) and tetragonal (T) phases, given 
as a fraction of a or c, with respect to ideal positions (i.e., for 
S x , with respect to the M x symmetry plane). 



Atom 


S z (C) 


8. (T) 


MT) 


Ba(l) 


-0.0279 


-0.0142 


-0.0277 


Om(l) 


-0.0140 


-0.0298 


-0.0126 


Ti(2) 


0.0092 


-0.0086 


0.0098 


Oi(2) 


0.0048 


-0.0297 


0.0059 


0„(2) 


0.0048 


-0.0240 


0.0045 


Ba(3) 


-0.0053 


-0.0149 


-0.0059 


Om(3) 


-0.0026 


-0.0280 


-0.0020 


Ti(4) 





-0.0034 





Oi(4) 





-0.0340 





0„(4) 





-0.0256 






Ti02 in the relaxed rutile structure. 

Now, we define the two chemical potentials /XTi0 2 an d 
MBaO in such a way that /iTi0 2 = corresponds to a sys- 
tem in contact with a reservoir of bulk crystalline Ti02, 
and similarly for ^BaO- Furthermore, if we insist that the 
system is always in equilibrium with a reservoir of bulk 
BaTiOa, then we have that 



fJ-BnO + MTi0 2 



-E f 



(2) 



Thus, only one of /iBaO and /iTi0 2 is an independent de- 
gree of freedom. We arbitrarily chose /iTio 2 as the inde- 
pendent one. Then /uxi0 2 can be allowed to vary over the 
range 



Ef < /i T i0 2 < 



(3) 



At /iTi0 2 = —Ef the system is in equilibrium with BaO 
and BaTiOs, and for lower values bulk crystallites of BaO 
can precipitate. Similarly, above [iTiOi = 0j bulk crys- 
tallites of Ti02 can form. 

Therefore, the grand thermodynamic potential per sur- 
face unit cell is given by 



F — y [Eslab - A r T i0 2 (MTiQ 2 + ^TiQ 2 



--^BaO^BaO + -E-Bao)] 



(4) 



(the factor of 1/2 appearing because the cell contains two 
surfaces), where -E s i a b is the energy of the relaxed slab in 
the tetragonal phase. For example, for the type-I slab, 
one has A Ti o 2 = 3 and A Ba o = 4. Eqs. (|2|) and (||) give 
F as a function of ^Ti0 2 over the range of Eq. (||). 



III. RESULTS AND DISCUSSIONS 

A. Structural relaxations 

First we determined the equilibrium atomic positions 
for our two types of slabs in the two phases. For the 



TABLE III. Atomic relaxations of the Ti-terminated sur- 
face (slab II) in the cubic (C) and tetragonal (T) phases, given 
as a fraction of a or c, with respect to ideal positions. 



TABLE IV. Calculated interlayer relaxation (/?) and rum- 
pling (77) for the surface layer of the relaxed slabs in the cubic 
(C) and tetragonal (T) phases (A) . 



Atom 



5 Z (C) 



5 X (T) 



*. (T) 



Slab 



/3(C) 



V(C) 



13 (T) 



»?(T) 



Ti(l) 


-0.0389 


0.0005 


-0.0331 


Oi(l) 


-0.0163 


-0.0499 


-0.0100 


On(l) 


-0.0163 


-0.0366 


-0.0071 


Ba(2) 


0.0131 


-0.0148 


0.0186 


Om(2) 


-0.0062 


-0.0292 


-0.0023 


Ti(3) 


-0.0075 


0.0019 


-0.0058 


Oi(3) 


-0.0035 


-0.0372 


-0.0022 


On (3) 


-0.0035 


-0.0278 


-0.0023 


Ba(4) 





-0.0111 





Om(4) 





-0.0276 






cubic surface, we started from the ideal structure and re- 
laxed. For the tetragonal surface, we obtained a starting 
guess by superposing the z-displacements from the re- 
laxed cubic surface with ^-displacements from the bulk 
tetragonal structure.^ The relaxed geometries are sum- 
marized in Tables H and [II. In these tables, the atoms 
are listed in the same order as shown in Fig. 1. (Coor- 
dinates are only listed for atoms in the top half of the 
slab, z > 0; the others are determined by the M z mirror 
symmetry.) By symmetry, there are no forces along x 
or y for the cubic surfaces, and no forces along y for the 
tetragonal surface. Also due to the crystal termination, 
the two O atoms associated with the Ti atom (Oi and 
On), are no longer equivalent in the tetragonal phase (I, 
II, and III indicate the O that is connected to Ti by a 
bond along x, y, and z, respectively). 

From Tables || and III, we can see that the largest 
relaxations are on the surface-layer atoms, as expected. 
They are especially important for the Ti atoms in the 
Ti-terminated slabs, in the direction perpendicular to 
the surface. This can plausibly be explained by noting 
that for bulk BaTiOa, the filled Ba levels lie well be- 
low the oxygen-2p valence-bands and do not hybridize 
strongly, so that the Ba atom is a relative spectator in 
the bondingEa In the tetragonal case, we can also see 
that the asymmetry between the two O atoms lying in a 
common surface plane is significant. Bond lengths on the 
surface change by less than 1.5% with respect to the bulk 
in the same phase (the latter yaJues are the experimental 
ones, taken from Mitsui et alE3). 

We computed the average displacements /3 and the 
rumpling 77 for the surface layers; the results are given in 
Table IV. To fix notation, let Sz(M) be the change in the 
surface-layer metal-atom z position relative to the ideal 
unrelaxed structure, and 5z(0) be the same for the sur- 
face oxygens (defined as [6z(0{) + <5z(On)]/2 for a Ti0 2 
layer). Then the surface relaxation parameter (3 is de- 
fined as [dz(M) + Sz(0)]/2, and the rumpling r\ is defined 
as [6z(0) — 5z(M)]/2. We find that the surface layers con- 
tract substantially inwards towards the bulk, with both 



Slab I 
Slab II 



-0.08 
-0.11 



0.03 
0.05 



-0.08 
-0.08 



0.03 
0.05 



the metai-and oxygen ions relaxing in the same direction. 

CohenE-D has computed the surface relaxations of an 
asymmetrically-terminated cubic slab (BaO on one sur- 
face and TiC>2 on the other) . Detailed quantitative agree- 
ment is probably not expected, because (i) only the 
surface-layer atoms were relaxed in Cohen's calculation, 
and (ii) the asymmetric termination introduces a small 
electric field which may have influenced the relaxations. 
Nevertheless, we do find qualitative agreement. Cohen 
finds that the Ba and O atoms relax inwards by 0.043 
and 0.033 lattice constants, respectively, on the Type-I 
surface; and the Ti and O atoms relax inwards by 0.048 
and 0.027, respectively, on the type-II surface. These 
can be compared with th e fir st two entries in the S Z (C) 
column of Tables || and III. It can be seen that Co- 
hen's values are 20-100% larger in magnitude, but of the 
same sign, as those that we calculate. Similarly, Cohen 
computes values of -0.15 A for the average surface layer 
relaxation, to be compared with the values given by us in 
the column /3(C) of Table IV. The rumpling computed 
by Cohen is also of the same sign, but different in detail, 
as that calculated by us. 

We are not aware of any experimental surface structure 
determination for BaTiC>3 with which we can compare 
our theory. However, we note that.a-LEED I-V study of 
the corresponding SrTi03 surfacestil comes to an oppo- 
site conclusion, suggesting an outward relaxation of the 
surface layer on the order of 0.1 A. While we have not 
carried out a parallel calculation on the SrTiOs surface, 
it appears unlikely that the mere replacement of Sr by 
the chemically similar Ba could be responsible for such a 
large qualitative discrepance- Thus, we suggest that the 
experimental interpretationEil should be reexamined. 



B. Influence of the surface upon ferroelectricity 

It is important to understand whether the presence 
of the surface has a strong effect upon the near-surface 
ferroelectricity. For example, is the FE order enhanced 
near t he sur face, or is it suppressed? As we shall see 
in Sec. IIIC , the energy scale of the surface relaxations 
is larger than the energy scale of the FE double-well 
potential. Thus, a strong effect is possible. To an- 
alyze whether it really occurs, we computed an aver- 
age FE distortion Sfe for each layer of the slab. We 
define <5fe = ^(Ba) — S x (Om) for a BaO plane, and 
^fe = <y a (Ti) - [4(00 + 4(0„)]/2 for a Ti0 2 plane. 



TABLE V. Calculated FE distortion <5 FE of the relaxed 
slabs, for each layer (units of lattice constant). Last line gives 
theoretical bulk values for reference. 

Slab I Slab II 

Layer <5 FE (BaO) <5 FE ( Ti 2 ) <5 FE (BaO) <5 FE (Ti0 2 ) 



1 


0.0131 






0.0345 


2 




0.0264 


0.0165 




3 


0.0157 






0.0434 


4 




0.0259 


0.0144 




Bulk 


0.0232 


0.0278 


0.0232 


0.0278 



Our calculated values for <5fe are given in Table [V|; the 
last row of the table gives the bulk values for reference. 
No clear pattern appears to emerge from these results, 
although we do note a moderate enhancement of the FE 
instability in the TiC>2 layers for the Ti02-terminated 
surface. 

The lack of a clear trend for the influence of surface 
effects upon the FE distortion can be understood, at least 
in part, by noting that the FE mode is only one of three 
zone-center modes having the same symmetryO The FE 
mode is distinguished as the one that is soft (ui 2 < 0) 
in the cubic structure, but there is no particular reason 
why the surface relaxation should couple more strongly 
to this mode than to the others. We have estimated how 
strongly the surface relaxations are related to each of 
the zone-center modes by the following procedure. We 
calculate the forces for a tetragonal surface slab in which 
the displacements in the x direction are those of the ideal 
bulk tetragonal structure, while the displacements in the 
z direction are taken from the relaxed cubic surface. The 
forces in the x direction are then projected onto each of 
the zone-center bulk modes polarized along x. That is, 
the force for each type of atom (e.g., Om) was summed 
over all such atoms in the slab, and the inner product 
was then taken between the resulting force vector and 
the bulk mode eigenvectors. 

We found that the FE mode accounts for only about 
31% and 26% of this force vector for the type-I and type- 
II slabs respectively. Thus, it seems that the distortions 
induced by the presence of the surface are mostly of non- 
FE character, helping to explain why the FE order is not 
as strongly affected as might have been guessed. 



C. Surface energies 

We turn now to a study of the surface energetics. Fol- 
lowing the approach outlined in Sec. ||, we calculated the 
grand thermodynamic potential F for our two types of 
surface as a function of TiC>2 chemical potential. The re- 
sults are shown in Fig. 2. In order to attain optimal can- 
celation of errors, £Tio 2 an d -^BaO were calculated within 
the LDA using the same pseudopotentials, and with the 
same 25 Ry energy cut-off. A similar k-point sampling 
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FIG. 2. Grand thermodynamic potential Fasa function 
of the chemical potential /iTi0 2 1 f° r the two types of surfaces, 
in the tetragonal phase. Dashed and solid lines correspond 
to type-I (BaO-terminated) and type-II (Ti02-terminated) 
surfaces respectively. 



as for the surfaces of BaTi03 was also employed. In 
this way, we obtained E{ = 3.23 eV for the formation en- 
ergy of BaTiC>3. This quantity fixes the range of physical 
values of //Ti0 2 > the left and right boundaries of Fig. 2 
correspond to a system in thermodynamic contact with 
bulk BaO and bulk TiC>2 respectively. It can be seen 
that both surfaces have a comparable range of thermo- 
dynamic stability, indicating that either type-I or type-II 
surfaces could be formed depending on whether growth 
occurs in Ba-rich or Ti-rich conditions. 

The average surface energy E sur { (i.e., the average of 
F for the two kind of surfaces) is independent of /iTio 2 ■ 
Therefore, this quantity is suitable for comparisons. For 
the cubic phase, we estimated E sur f for the (001) surfaces 
to be 1.241 eV per surface unit cell (1265 erg/cm 2 ); and 
for the tetragonal phase, it was estimated to be 1.237 
eV per surface unit cell (1260 erg/cm 2 ). The value of the 
average E SUI f calculated in Rcf. ttj\ for the symmetrically- 
terminated cubic (001) surfaces is 920 erg/cm 2 . As 
pointed out in that paper, the large value of E SUT f may 
help explain why BaTiC>3 does not cleave easily, but frac- 
tures instead. 

In order to compute the surface relaxation energy 
^rciax, we computed the average surface energy E unie \ for 
the unrelaxed cubic slabs (i.e., atoms in the ideal cubic 
perovskite positions), using the same k-point sampling as 
for the relaxed systems. We obtained E unrc \ = 1.358 eV. 
Thus, the relaxations account for 0.127 eV of the surface 
energy per surface unit cell (or about 130 erg/cm 2 ). 

Note that -E r dax is many times larger than the bulk 
ferroelectric well depth, estimated to be of the order of 
0.03 eV. This would indicate that the surface is capa- 
ble of acting as a str ong p erturbation on the FE order. 
As explained in Sec. IIIB, however, the actual effect is 
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TABLE VI. Calculated band gaps for relaxed cubic (C) 
and tetragonal (T) surface slabs (eV). 
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FIG. 3. Calculated band structures for BaTi03 in the 
cubic phase, (a) Surface-projected bulk band structure, (b) 
BaO-terminated surface (slab I), (c) Ti02-terminated surface 
(slab II) . The zero of energy corresponds to the bulk valence 
band maximum. Only the lowest few conduction bands are 
shown. 



more modest than one would guess based on energetic 
considerations alone. 



D. Surface band structure 

We have carried out LDA calculations of the surface 
electronic structure for our various surface slabs. While 
the LDA is well known to be quantitatively unreliable 
as regards excitation properties such as band gaps, we 
believe that the results presented here are nevertheless 
likely to be qualitatively correct. The bulk band gap 
in our calculation is 1.8 eV, to be compared with the 
experimental value of 3.2 eV; this level of disagreement 
is typical for the LDA. 

Fig. 3 shows the calculated LDA band structure for 
cubic bulk BaTiC>3 projected onto the surface BZ, and 
the surface band structures for the Ba- and Ti-terminatcd 
relaxed surfaces in the cubic phase. (The zero of energy 
for each surface slab was established by aligning the Ba 
or Ti semicore s states in the interior layers of the slab 
with those of the bulk.) Plots for the tetragonal surfaces 
would look similar, except that the tendency for states to 
intrude into the gap is stronger for the cubic case. This 



Slab 



C 



T 



Slab I 
Slab II 
Bulk 



1.80 
0.84 
1.79 



2.01 
1.18 
1.80 



can be seen in Table VI , where we list the calculated band 
gaps for both cubic and tetragonal slabs. We therefore 
tend to focus on the cubic surfaces, where it is easier to 
identify and characterize the surface states. 

First, we can see that on the Ba-terminated (type-I) 
surface, the gap is not reduced and there are no deep gap 
states. On the Ti-terminated (type-II) surface, however, 
the gap is reduced, especially for the cubic case. As can 
be seen from Fig. 3(c), there is a tendency for valence- 
band states to intrude upwards into the lower part of the 
band gap for this surface, especially near the M point 
of the surface BZ. (Qualitatively similar results can be 
seen in Fig. 3 of Ref. 17.) However, the conduction band 
does not change much with respect to the bulk. More- 
over, we do not see signs of any true "deep-gap" states 
lying near the center of the gap. As noted earlier, the 
existence of such deep-gap states remains controversial 
experimentally. Our work suggests that if gap states do 
exist in connection with non-defective (001) surfaces, it 
is likely that they would be found in the lower part of the 
band gap, and that this would be indicative of exposed 
Ti02 (as opposed to BaO) surface planes. 

Figure 4 illustrates the character of the valence-band 
state at the M point that is intruding into the lower part 
of the gap. The total charge density is also shown for 
reference. It can be seen that this state is composed of O 
2p lone-pair orbitals lying in the surface plane. Further 
inspection shows that the special feature of this state is 
that the wavefunction has four nodal planes [(100), (110), 
(010), and (110)] intersecting at the Ti sites. This pre- 
cludes the presence of any Ti 3d character (in fact, any 
Ti character of angular momentum I < 4). In the bulk, 
the oxygen 2p orbitals are all hybridized with Ti 3d or- 
bitals to some extent, and the level repulsion associated 
with this hybridization pushes the energy location of the 
valence-band states downward in energy. Thus, the en- 
ergy of the unhybridized O 2p lone-pair surface state at 
the M point is left intruding into the lower part of the 
gap. 

This insight makes it possible to understand other fea- 
tures of the surface band structures as well. For exam- 
ple, on the BaO-terminated surface, every surface oxy- 
gen atom is directly above a Ti atom, and is strongly 
hybridized to it. Thus, there is no such tendency for the 
formation of surface states in this case. Returning to the 
Ti02-tcrminated surface, there seems to be a weaker ten- 
dency for the intrusion of a valence-band derived surface 
state at T. This state turns out to have a single nodal 




FIG. 4. Charge-density plots for relaxed Ti-terminated 
surface (slab II) in the cubic phase. Panels (a-b) show cuts 
on a vertical x-z BaO plane; panels (c-d) show cuts on the x-y 
surface Ti02 plane. Total charge density, (a) and (c); charge 
density of highest occupied state at M point of the surface 
BZ, (b) and (d). 



plane passing through the Ti site, so that hybridization 
with Ti 3d orbitals is only weakly allowed (by the break- 
ing of M z mirror symmetry across the surface plane) . 

We expect that these results would remain qualita- 
tively valid for other II- VI perovskites such as SrTiC>3 
or PbZrC>3. For I-V perovskites such as KNbC>3 and 
LiTaC>3, however, the (001) surface is non-stoichiometric, 
and the surface electronic structure would be expected to 
be quite different. 



IV. SUMMARY 

In summary, we have carried out LDA density- 
functional calculations of BaO- and Ti02-terminated 
(001) surfaces for cubic and tetragonal phases of BaTiC>3. 
By minimizing the forces on the ions, we obtained the re- 
laxed ionic structures. As would have been expected from 
the bulk electronic levels of BaTiC>3, the most impor- 
tant relaxations occur for the Ti02-terminated surfaces. 
There appears to be a modest tendency for the surface 
relaxations to enhance the FE distortion on that surface, 
although the situation is complicated by the fact that the 
relaxations excite modes other than the soft zone-center 
one. 

The free energies for the different surfaces were cal- 
culated as a function of the TiC>2 chemical potential. 



In particular, the average surface energy was found to 
be about 1260 erg/cm 2 . The surface relaxation energies 
were found to be around 10% of the total surface energy. 
In accord with previous theoretical reports, no mid-gap 
surface levels are found. But for the TiC>2 terminated 
surfaces, there is a substantial reduction of the bulk gap, 
especially for the cubic-phase surface. This reduction 
results from the intrusion of states of valence-band char- 
acter into the lower part of the band gap, especially near 
the M point of the surface BZ. 
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